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摘 要 : 作物 系数 (K.) 对 提高 实际 蒸 散 量 的 估算 精度 及 区 域 水 资源 调控 具有 重要 意义 。 通 过 Landsat 8 数据 计算 


= 
T 


归 一 化 植被 指数 (NDVI) .土壤 调整 植被 指数 (SAVI) 和 简单 比值 植被 指数 (SR) ,结合 2017 年 地 面 实测 土壤 含水 量 
(SM) 和 叶 面 积 指 数 (LADT) ,通过 逐步 回归 分 析 建 立 科 尔 学 


也 区 草 多 及 沙丘 试验 区 的 作物 系数 估算 模型 ,利用 2018 


年 模拟 作物 系数 值 与 FAO 56 P-M 模 型 计算 所 得 潜在 蒸 散 量 (ETu) 相 乘 ,得 到 实际 蒸 散 量 (ET,) 的 估算 值 ,与 涡 度 相 


关系 统 所 测 实际 蒸 散 量 进 行 对 比 验 证 。 结 果 表 明 :(1) 草 


多 和 沙丘 试验 区 的 作物 系数 在 生长 季 内 的 变化 趋势 与 植 


被 指数 及 SM 具有 一 致 性 ,说 明了 基于 这 些 指标 建立 玉 佑 算 模型 的 可 行 性 ;(2) 相关 性 分 析 中 , 草 多 试验 区 的 玉 与 
SR 相关 性 不 显著 (P>0.05 ) ,沙丘 试验 区 玉 与 SR 相关 性 较 低 (0.46) , 故 剔 除 该 因子 ;(3 ) 逐步 回归 分 析 中 进一步 剔除 
不 显著 因子 后 建立 草 旬 及 沙丘 作物 系数 估算 模型 ,模型 的 均 方 根 误差 修正 的 决定 系数 分 别 为 0.06 .0.84 和 0.12、 


0.71; (4) 经 涡 度 相关 数据 验证 ,基于 K. 佑 算 模 型 计算 的 87, 在 草 僻 和 沙丘 试验 区 均 取 得 较 好 的 模拟 效果 。 
遥感 ; 作物 系数 ; 植被 指数 ， 涡 度 
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地 表 蒸 散发 (Evapotranspiration, 有 ET) 由 土壤 蒸 
发 和 植被 蒸腾 组 成 ", 是 连接 生态 与 水 文 过 程 的 重要 
纽带 ,同时 也 是 地 表 水 热平衡 的 重要 组 成 部 分 。 
据 统计 ,全 球 约 60% 以 上 的 降水 以 蒸 散 的 形式 耗 散 ” ， 
而 在 干旱 半 干 旱地 区 ,以 蒸 散 形式 返回 到 大 气 系统 
中 的 水 分 甚至 达到 80% 以 上 ,因此 ,准确 的 合算 地 
表 蒸 散发 量 , 对 维持 干旱 半 干 旱地 区 生态 系统 平 
衡水 资源 合理 开发 利用 及 区 域 小 气候 调节 等 具有 
重要 意义 "。 

常用 的 燕 散 发 估算 方法 有 : 波 文 比 能 量 平衡 
法 ARB WOK S EAA PAO 56 作物 系数 法 等 ， 
其 中 ,FAO 56 作 物 系 数 法 是 1998 年 由 联合 国 粮农 
组 织 (Food and Agricultural Organization, FAO ) 基于 
Penman-Monteith 公式 ,利用 参考 作物 燕 散 量 (Refer- 
ence Evapotranspiration, ET)) 以 及 作物 系数 到 计算 
农田 1 的 方法 , FAO 56 Penman-Monteith 公式 将 
实际 作物 冠 层 阻力 和 空气 动力 阻力 与 参考 作物 的 
差异 均 包 含 在 作物 系数 到 中 ,能够 贴切 .简单 的 反 
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映 控制 葵 散 的 物理 和 生理 因素 ,对 7 的 模拟 精度 
高 .操作 简单 .实用 性 强 ,在 世界 范 于 内 都 得 到 广泛 
应 用 四。Taylor 等 "利用 FAO 56 作 物 系数 法 估算 了 
3 个 柑橘 园 的 水 分 利用 情况 ,结果 表明 :该 模型 具有 
一 定 适 用 性 ; Anderson 等 … 通 过 有 7 与 双 作 物 系 数 
相 乘 得 到 了 农业 灌区 的 蒸 散 量 ,并 用 涡 度 相关 进行 
验证 ,结果 显示 模拟 效果 较 好 ; 李 新 等 "通过 单 、 双 
作物 系数 法 计算 砾石 覆盖 条 件 下 小 麦 蒸 散 量 ,发 现 
模拟 值 与 实测 值 拟 合 良好 ;于 文 颖 等 ”研究 得 出 改 
进 后 FAO-56 模 型 的 双 作 物 系 数 法 可 以 有 效 的 估算 
辽河 三 角 洲 湿地 芦苇 群落 的 蒸 散 量 。 

由 于 作物 系数 具有 一 定 的 经 验 性 ,并 不 适用 于 
所 有 地 区 ,直接 使 用 推荐 的 作物 系数 估算 蒸 散 量 的 
误差 较 大 ,作物 系数 受 植被 生长 状况 、 气 候 及 水 分 
条 件 等 诸多 因素 的 影响 ,必须 结合 当地 实际 情况 进 
行 校正 下 。 但 这 些 因素 所 造成 的 影响 很 难 及 时 对 
作物 系数 曲线 作出 调整 ,而 植被 指数 能 够 实时 快速 
的 反映 作物 的 生长 状况 ,这 使 得 遥感 技术 在 玉 的 佑 
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算 中 得 以 广泛 应 用 。Duchemin 等 分析 了 通过 
Landsat 数据 获得 的 NDVI 来 间接 估计 作物 系数 的 可 
能 性 ; Glenn 4" AJ FH MODIS 卫星 遥感 影像 获取 的 
叶 面 积 指数 (LAT) . 归 一 化 植被 指数 (NDVI) .土壤 调 
整 植被 指数 (SAVJ) 来 估算 到, 可 以 很 好 的 模拟 预测 
从 沙漠 ,北极 昔 原 到 热带 雨林 等 多 种 生物 群 类 型 的 
ET; EESE X T 3b ifii SC Al Landsat 8 卫星 遥感 
反 演 的 叶 面 积 指数 ,建立 了 作物 系数 与 叶 面 积 指数 
的 经 验 对 数 模型 ,结果 显示 在 低 植被 覆盖 下 优化 后 
的 模型 精度 较 高 。 

科尔沁 沙 地 是 我 国 典 型 的 干旱 oe E MUSS 
错过 度 带 ,沙丘 和 草 甸 相间 分 布 ”2 , ET, AN DUCK 
文 循环 中 占 主导 地 位 ,同时 也 是 干旱 区 水 文 研 究 
的 重要 环节 之 一 。 虽 然 基于 植被 指数 的 及 估算 模 
型 取得 较 好 的 效果 ,但 作物 系数 不 仅仅 受 植被 生理 


加 作物 系数 模拟 及 蒸 散 发 估算 


1 研究 区 概况 


研究 区 位 于 内 蒙古 自治 区 通辽 市 科尔沁 左 恤 
后 旗 阿 十 拉 生 态 水 文 试验 站 (图 1) ,面积 约 55 
km”。 该 地 区 多 年 平均 降雨 量 389 mm, 主要 集中 在 
6 一 9 月 ;多 年 平均 气温 约 为 6.6 %C, 多 年 平均 相对 湿 
度 为 55.8%。 地 势 南北 高 .中 间 低 ,由 北向 南 形成 了 
典型 的 沙丘 -过 渡 带 - 草 甸 地 -湖泊 - 草 旬 地 -沙丘 的 
地 貌 形态 。 人 研究 区 属 温带 半 干 旱 大 陆 性 季风 气候 
区 , EK PR A RA 22 D TRES (Artemisia haloden- 
dron) ,/\\ HH $8 35 JL (Caragana microphylla) | Kf E5 
(Artemisia sieversiana) , i Wl (Salix gordejevii) Fil 1% TU 
(Populus )^$ , 

草 多 试验 区 (C4) 位 于 研究 区 中 间 地 势 较 低 区 
域 ,紧邻 湖泊 , 样 地 地 势 平 坦 , 地 下 水 埋 深 浅 , 土 壤 


生态 特征 影响 ,还 与 土壤 水 分 有 关 。 综 上 所 述 ,本 
文 利 用 Landsat 8 卫星 数据 ,提取 3 种 植被 指数 : 归 一 
化 植被 指数 (NDVI) .土壤 调整 植被 指数 (SAVI) 和 
简单 比值 植被 指数 (SR)22 ,结合 地 面 实测 表层 土 
培 含 水 量 (SM) 及 叶 面 积 指数 (LAI) ,分 别 建 立 草 负 
及 沙丘 试验 区 玉 佑 算 模 型 ,并 用 涡 度 数据 对 模型 的 
可 行 性 及 适用 性 进行 评价 ,为 作物 系数 在 区 域 尺度 
上 的 应 用 及 和 干旱 半 和 干旱 地 区 葵 散 量 的 准确 估算 提 
供 科学 依据 。 
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图 1 研究 区 及 试验 点 位 置 示意 图 


长 期 处 于 湿润 状态 ,土质 为 沙壤土 ,主要 植被 有 芦苇 
(Phragmites communis ) 4 (Leymus chinensis)、 灯 让 
38 (Juncus effusus) 等 ;沙丘 试验 区 (A4) 位 于 研究 区 北 
部 的 沙丘 地 带 , 植 被 覆盖 度 较 低 ,周围 生长 的 天 然 植 
BAF GIS Wx (Agriophyllum squarrosum) 等 。 


2 数据 与 方法 


2.1 数据 采集 与 处 理 
2.1.1 气象 土壤 及 涡 度 数据 草 旬 和 沙丘 试验 区 
各 安装 有 一 套 开路 涡 度 相关 系统 和 和 气象 观测 塔 , 涡 


流动 沙丘 he EEA a NH 
。 草 向 湿 地 试验 点 号 典型 区 边界 


wae Lx 
NE mm 农田 


Fig. 1 Locations and sampling sites of the study area 
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度 相 关系 统 包 括 开路 红外 气体 分 析 仪 .三维 超声 风 
速 仪 以 及 四 分 量 净 辐射 仪 .光合 有 效 辐 射 仪 ,测定 
LEH fil CO. 38 E , 以 及 净 辐 射 . 光 合 有 效 辐射 等 要 
素 ,采样 频率 为 30 Hz; 波 文 比 气象 土壤 环境 监测 系 
统 包 括 温 湿度 测量 仪 .风速 风向 测定 仪 和 翻斗 式 雨 
He fa) ,主要 观测 的 气象 数据 包括 空气 温 湿 度 .风速 、 
风向 、 降 雨量 等 ;土壤 环境 要 素 可 测定 土壤 热 通 量 、 
土壤 温度 .土壤 含水 量 等 ( 表 1)。 

利用 美国 Li-COR 公司 开发 的 Eddy Pro 软件 ,对 
涡 度 相关 系统 采集 的 数据 进行 处 理 , 包 括 二 次 坐标 
旋转 校正 ,频率 损失 修正 、 趋 势 校 正 .、WPL 校 正 等 一 
系列 修正 。 通 量 数据 质量 控制 主要 包括 降水 期 间 
数据 的 剔除 .夜间 数据 的 筛选 与 剔除 . 闪 值 检验 与 
野 点 剔除 和 大 气 层 结 状 况 下 的 数据 剔除 。 为 保证 
数据 的 完整 性 ,小 于 3 的 数据 采用 线性 内 插 法 进 
行 插 补 ; 较 长 时 间 的 缺失 数据 则 采用 平均 每 天 变异 
法 (MDV)* 进 行 插 补 ,进而 得 到 较为 精准 的 通 量 数 
据 。 草 旬 和 沙丘 试验 区 的 能 量 平衡 闭合 率 EBR 分 
别 为 70% 和 76%, 通 量 数据 可 信和 度 高 ,能够 代表 观测 
区 域 通 量 大 小 。 
2.1.2 植被 数据 在 涡 度 观测 源 区 内 ,利用 Li-2200 
叶 面 积 指数 仪 测定 LAI。 草 旬 试 验 区 植被 长 势 密集 
且 均 一 ,以 气象 站 为 中 心 共 布 设 9 个 草本 样 方 (1 mx 
lm ) 进 行 调查 ;沙丘 区 植被 覆盖 度 较 低 , 共 布 设 3 个 
30 mx30 m 的 大 样 方 ,每 个 大 样 方 中 选 取 对 角 线 及 
中 心 的 5 个 的 半 灌 木 样 方 (3 mx3 m)。 以 上 样 方 的 
选取 均 可 代表 涡 度 源 区 内 植被 指标 的 平均 水 平 ,以 
Landsat-8 重 访 周 期 为 调查 周期 ,每 月 2 次 ,同时 记 
录 样 方 内 植被 株 高 等 指标 。 
2.1.3 RAYE 采用 的 遥感 数据 为 Landsat8 OLI 
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数据 。 在 数据 选择 方面 ,为 准确 计算 植被 指数 , 根 
据 研 究 区 已 有 卫星 数据 的 影像 质量 、 成 像 时 间 等 信 
县 ,选择 2017 一 2018 年 生长 季 4 一 10 月 的 19 期 无 云 
或 少 云 ,清晰 度 高 的 遥感 影像 ,影像 数据 取 自 美国 
地 质 调查 局 (USGS) 网 站 (http:/earthexplorer.usgs. 
gov/) 以 及 地 理 空间 数据 云 网 站 (http://www.gscloud. 
cn/)。 采 用 ENVI5.3 软件 对 获取 的 遥感 影像 进行 几 
何 校正 辐射 定 标 直方 图 均衡 化 .研究 区 影像 裁剪 
等 预 处 理 , 以 提高 植被 指数 计算 精度 。 

2.2 研究 方法 

2.2.1 作物 系数 计算 方法 “作物 系数 (K) 是 实际 蒸 
散 量 (E7,) 与 参考 作物 蒸 散 量 (ET7) 之 比 ,综合 反 映 
了 作物 种 类 .生长 状况 .土壤 水 肥 条 件 等 因素 对 作 
39 Fe Bi ny se n °°) ,利用 FAO 56 Penman-Monteith 模 
型 输入 气象 数据 计算 E77,, 人 然后 根据 涡 度 塔 测 得 的 
潜 热 通 量 数据 与 &7, 的 比值 来 计算 K.: 


ET. 
K= RT, "y 
900 
0.408A(R, - G) +y ue, e, 
- ( ) T +273 ( ) (2) 


d A+y(1+0.34u,) 
SU KCN TE BM ET, Ab ZR HG (mm d') ; 
ET) HBS EIZE BCR (mm d); ACA AIK AUR TR 
FE HAR ES (KPa - 9C) ;及 为 冠 层 表面 净 辐 射 (MJ. 
md ') ;7 为 平均 气温 (°C);G6 为 土壤 热 通 量 (MJ- 
med!) ;也 为 2 m 高 处 风速 (ms ');e, 为 饱和 水 汽 压 
(KPa) ;e, 为 实际 水 汽 压 (KPa) ;7 为 湿度 计 和 常数 
(KPa-9C7), 
222 植被 指数 计算 方法 “植被 指数 是 不 同 遥 感光 
谱 波 段 间 的 线性 或 非 线性 组 合 , 主要 反映 在 可 见 光 、 


表 1 地 面 监测 基本 信息 
Tab. 1 Generalinformation of the ground stations 
类 型 观测 项 仪器 型 号 eee ee ie 采样 频率 / 
草包 沙丘 (min: 1X7) 
SUE 开路 红外 气体 分 析 仪 Li-7500,USA 2.7 4.95 30 
相关 三 维 超声 风速 CSAT-3 ,USA 
系统 辐射 分 量 NR-LITE,USA 3.5 3.5 
空气 温 湿度 HMP45C 1.5 .2.4 2.5.10 10 
气象 风速 风向 034B 2.3.5.10 2.10 
~ 降雨 量 TE525MM 0.7 0.7 
BR 土壤 温度 ,水 分 HydraProbeCS655 0.05 .0.1 .0.2 0.1.0.2,0.4,0.6 
土壤 热 通 量 HFPO1sc 0.1 0.1,0.2,0.4 
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近 红 外 波段 反射 与 土壤 背景 之 间 的 差异 ,可 定量 说 
明 植被 的 生长 状况 ”。 为 了 找到 与 Xk. 相关 性 最 好 的 
植被 指数 ,本 文选 取 归 一 化 植被 指数 (NDVI) .土壤 
调整 植被 指数 (SAVI1) 及 简单 比值 植被 指数 (SR)( 表 
2) , 草 负 和 沙丘 试验 区 的 各 植被 指数 通过 统计 涡 度 
源 区 内 所 有 像 元 的 均值 所 得 。 


表 2 植被 指数 计算 


Tab. 2 Vegetation indices used in this study 


名 称 公式 
r # Prin ~ Prev 
ES H NDVI = Pv Prev 
A INE TRUM ; Pxn + Pre 
1+L "e 
土壤 调整 植被 指数 (SAVD — savis LED o Puo) 
L * (Dun *puy) 
简单 比值 植被 指数 (SR) sR- Pm 


注 :2DNIR 为 近 红 外 波段 ;pn 为 红 波段 3 LIX 0.5. 


22.3 模型 一 致 性 评价 ”利用 2017 年 生长 季 气 象 、 
土壤 及 植被 数据 建立 K. 估 算 模 型 ,将 2018 4E Land- 
sat 8 数据 佑 算 的 植被 指数 、 地 面 实 测 叶 面积 指数 
(LAT) 及 表层 土壤 含水 量 (SM) 数 据 输入 估算 模型 来 
模拟 K., 并 与 FAO 56 P-M 模型 计算 所 得 BT, 相 乘 得 
到 7 的 模拟 值 ,通过 与 涡 度 相关 实测 实际 蒸 散 量 
进行 对 比 来 判别 玉 佑 算 模型 的 可 行 性 。 选 取 相 对 
均 方差 (RMSE) 一致 性 指数 (d) .平均 偏差 (bias) 和 


1P—s—NDVI —*x—SAVI 一 上 SR —e—LAI —e—K, ---4—- SM Ef, 
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4] E] RR ZEA A: 


Nash-Sutcliffe 效率 指数 (E)4 个 统计 量 ”对 模型 的 
精度 进行 评价 ,计算 公式 为 : 


RMSE = ade - oy (3) 


rer 
Y(P.-0) 
d21-|4———— (4) 
Y (P; 0...] *|0;* Onl) 
bias = 150 =P) (5) 
> 
E s (6) 
2.(0,- 0 ) 


3X HB iN D ULTRI CC; P: A SE E ZS Be IE I [EE 
(mm: d~’) ; 0,29 SEE ZI E E HY SC WU EL (mm * d) ; Pave 
为 时 间 尺 度 内 模拟 值 的 平均 值 (mm…d 7) ;0 为 时 间 
尺度 内 实测 值 的 平均 值 (mmd"')。 


3 结果 与 分 析 


3.1 相关 因子 动态 变化 分 析 
为 确定 各 因子 作为 自 变 量 的 可 行 性 ,首先 分 析 
了 2017 年 草 甸 及 沙丘 试验 区 的 各 植被 指数 (NDVI、 


Xx 


SAVI SR LAI) KEREKE (SM) KRAE 


ET, 


s : 
05-01 05-16 05-31 06-15 06-30 07-15 


07-30 | 08-14 08-29 09-13 09-28 


日 期 /月 -日 


注 :P 表 示 降 水 量 。 下 同 。 
图 2 2017 年 草 务 试验 区 相关 参数 变化 趋势 


Fig. 2 The variation trend of relevant parameters in the meadow test area in 2017 
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CET.) HERRE CET) TE I CK ) ERNE 
长 季 的 动态 变化 。 

由 图 2 可 知 ,2017 年 草 旬 试 验 区 生长 季 ET, ON 
528 mm, HJ 有 7 为 3.1$ mm ,最 大 值 出 现在 7 月 , 达 
6.47 mm, 最 小 值 在 5 月 ,为 0.99 mm, ERA MET) 
为 825 mm; NDVI,SAVI,SR 及 LAI 的 变化 趋势 大 致 
相同 ,生长 初期 变化 较 小 (5 月 初 ) ,生长 旺盛 期 随 着 
植被 的 生长 ,植被 覆盖 度 变 大 ,各 植被 指数 快速 增 
长 , 因 5 一 7 月 降水 较 少 ,6 月 底 植被 指数 有 所 下 降 ,8 
月 雨水 充沛 ,植被 生长 旺盛 ,但 由 于 在 8 月 16 日 进 
行 了 人 工 浊 割 牧草 ,各 植被 指数 急剧 减 小 ,之 后 略 
有 回升 ,9 月 中 旬 以 后 ,植被 进入 生长 末期 ,植被 指 
数 旺 下降 趋势 ;K. 与 植被 指数 有 相似 的 变化 规律 , 生 
长 旺 戌 期 高 于 生长 初期 和 末期 ,相关 学 者 在 以 往 的 
研究 中 也 证 实 K. 与 植被 指数 变化 规律 相 吻 合 *”; 
SM 在 生长 初期 有 明显 的 上 升 趋势 ,生长 中 期 维持 较 
高 水 平 ,而 在 生长 末期 由 于 植被 覆盖 度 降 低 , 土壤 
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而 SR 在 整个 生长 季 变 化 规律 不 显著 ,这 是 由 于 SR 
主要 用 于 高 密度 植被 区 域 的 动态 监测 3;K. 亦 和 植 
被 指数 呈现 相似 变化 趋势 , 即 先 增加 后 减少 ;SM E 
化 波动 较 小 ,生长 季 均 不 超过 10% ,高 的 SM 总 伴随 
着 降水 事件 的 发 生 。 

对 比 草 个 和 沙丘 试验 区 可 以 看 出 ,降水 不 仅 对 
土壤 水 分 产生 影响 ,同时 也 会 引起 K. 和 植被 指数 的 
变化 ,Xu 等 也 得 到 类 似 的 研究 结果 。 降 水 对 植被 
指数 的 影响 主要 表现 在 持续 干旱 后 植被 指数 有 减 
小 的 趋势 , 彭 文 甫 等 "研究 表明 ,植被 指数 与 降水 
量 呈 正 相 关 关 系 ; 草 旬 试 验 区 临近 湖泊 ,地 下 水 埋 
深 较 浅 ,土壤 水 分 较 大 ,K. 变 化 主要 受 植 被 生长 时 期 
的 影响 ,而 沙丘 试验 区 的 K. 对 降水 更 为 敏感 ,主要 
表现 在 K. 旦 波动 性 增长 , 雨 后 K 有 明显 的 上 升 趋 
势 ,而 连续 干旱 后 ,土壤 和 植被 都 处 于 水 分 气 缺 状 
态 , 蔡 散 量 减 小 进而 导致 作物 系数 呈 下 降 趋 势 , 生 


蒸发 加 强 ,SM 出 现 明 显 的 回落 。 

由 图 3 可 知 ,沙丘 试验 区 薰 散 量 大 致 呈 钟 形 曲 
线 分 布 , 生 长 季 内 总 1 为 247 mm, 生 长 初期 .生长 
旺盛 期 及 生长 末期 的 ET 分 别 为 30 mm ,171 mm 及 
46 mm ,分 别 占 总 本 ,的 12.15% .69.26% 和 18.61% , 
生长 季 内 总 ET 为 762 mm,5 月 ET 总 量 最 大 ,10 月 
最 小 ;NDVI、SAVI 及 LAI 变 化 趋势 一 致 ,各 值 整体 小 
于 草 多 试验 区 ,最 大 值 出 现在 8 月 ,最 小 值 在 4 月 ， 


长 过 程 和 降水 过 程 同 时 控制 着 生长 期 作物 系数 的 
变化 趋势 ,这 与 张强 等 3 的 研究 结果 一 致 。 
32 KK. 估 算 模 型 的 建立 

为 进一步 确定 .估算 模型 的 自 变 量 ,采用 相关 
性 分 析 筛 选 输入 因子 ( 表 3)。 草 甸 试 验 区 K. 与 ND- 
VI、SAVI、LAI 及 SM 的 相关 系数 分 别 为 0.74、0.69、 
0.78 及 0.80, 相 关 性 较 高 ,而 K. 与 SR 相关 性 不 显著 
(P>0.05) ;沙丘 试验 区 KK. 与 各 变量 相关 性 均 达 到 显 
著 水 平 ,K. 与 SR 的 相关 系数 为 0.46, 与 其 他 变量 相 


0 [人 + 
04-16 05-06 


07-25 


CF AR 


à 人 人 于 
09-03 09-23 10-13 


日 期 /月 -日 


图 3 2017 年 草 多 试验 区 植被 .土壤 及 蒸 散发 相关 参数 变化 趋势 


Fig. 3 The variation trend of vegetation, soil and evapotranspiration related parameters in meadow test area in 2017 


202012.00037v1 


chinaXiv 


李 BS CET DUE GE 


ChinaX iv EAT 


JEW RSA BB 


X3 作物 系数 K. 与 各 指标 间 相 关系 数 


Tab.3 Correlation coefficients between K. and each factors 


on 相关 系数 
站 点 

NDVI vs K. SAVI vs K. SR vs K. LAI vs K. SM vs K. 
草 人 多 0.74 0.69* 0.53 0.78" 0.80" 
沙丘 0.77” 0.74 0.46 0.85” 0.70 


注 :* 表 示 在 P<0.05 KF E T RE P<0.01 水 平 上 显著 。 


比 相关 性 较 弱 ,因此 ,本 文选 择 与 K. 相 关 性 最 高 的 
NDVI,SAVI,LAI & SM 这 4 个 变量 进行 估算 模型 的 

在 上 述 分 析 的 基础 上 ,利用 2017 年 生长 季 数 
据 , 将 4 个 变量 与 到 建立 逐步 回归 方程 ,逐步 回归 分 
析 的 基本 原理 是 将 变量 逐个 引入 到 模型 当中 ,对 每 
一 个 自 变 量 进行 假设 检验 , 当 某 一 个 自 变 量 对 K. 影 
响 不 显著 时 ,将 其 剔除 ,此 过 程 反复 进行 ,直至 既 无 
显著 的 变量 选 和 人 回归 方程 ,也 无 不 显著 自 变 量 从 回 
归 方 程 中 剔除 为 止 ,筛选 出 显著 影响 的 因子 作为 自 
变量 建立 最 优 归 回 方程 ”。 在 逐步 回归 分 析 过 程 
中 , 草 多 试验 区 将 不 显著 因子 (P>0.05)SAVI LAI ži 
除 ,沙丘 试验 区 剔除 NDVI 和 SAVI, 回 归 方 程 见 表 4。 

由 表 4 可 知 , 草 甸 试 验 区 利用 单一 的 NDVI 
(casel) fll NDVI, SM 组合 (case2) 均 取得 恨 好 的 效 
果 ,RMSE 分别 为 0.12 fH 0.06, XW K. 14] ffi Jj 0.65 
(0.41 ^ 0.89) , casel $X WW K. 35] (E Jy 0.64 (0.46 ~ 
0.81) , case2 模拟 均值 为 0.65(0.43 ~ 0.88) ; 沙丘 试 
Js EX K.-LAI(casel ) fl K.-LAI SM (case2) f] RMSE 分 
别 为 0.14 810.12, XM KIEA 0.36 (0.11 ~ 0.62) , 
casel ,case2 30135] (8.4) 91 Jy 0.35 (0.16 ~ 0.56) .0.36 
(0.17 40.61) 。 草 甸 和 沙丘 试验 区 的 case2 模拟 值 
更 接近 于 实测 值 , 另 外, 随 着 变量 的 加 入 , Roa? AT VA 
更 直观 的 看 出 ,模型 的 效果 不 断 提 高 , 草 多 试验 区 
HY) Roa? FH 0.52 提高 到 0.84, 沙丘 试验 区 由 0.57 提高 
到 0.71,SM 的 加 入 使 得 模型 效果 提升 显著 ,这 与 张 
瑜 等 ”的 研究 结果 一 致 。 

作物 蒸 散 是 土壤 和 植被 水 分 扩散 并 输送 到 大 


气 的 过 程 , 其 影响 因素 不 仅 包括 温度 .饱和 水 汽 压 
差 .风速 等 局 地 气候 条 件 ,还 受 限于 土壤 和 植被 等 
环境 因素 *”。NDVI 和 LAI 反 映 了 植被 的 物候 特 
征 ,包括 植被 的 密度 以 及 作物 生长 时 期 和 作物 高 度 
的 变化 等 ,可 以 有 效 地 描述 K. 在 生长 季 的 变化 趋 
势 。 土 壤 水 分 是 土壤 蒸发 和 植被 蒸腾 最 直接 的 水 
分 来 源 ,土壤 水 分 状况 (尤其 是 表层 ) 将 会 直接 影响 
地 - 气 界 面 水 气压 差 和 植被 水 分 散失 的 主要 通道 
(气孔 ) 的 开 度 ,是 反映 植物 水 分 胁迫 和 土壤 表面 水 
分 蒸发 的 指标 。 以 往 的 研究 多 基于 一 种 植被 指数 
建立 简单 线性 回归 方程 ,如 Gontia 等 中 建立 了 KK.- 
NDVI 和 K.-SAVI 模 型 对 印度 地 区 干渠 灌溉 指挥 部 
的 小 麦 作物 的 实际 作物 燕 散 量 进行 了 舍 算 ;Lei 等 
提取 了 NDVI、EVI 等 5 种 植被 指数 ,分 别 与 K. 建 立 
线性 方程 ;Zolfagharnejad 等 中 采用 GIS 软件 提取 
NDVI、DVI、SAVI 等 5 种 植被 指数 建立 冬小麦 K. 售 
算 模型 ,结果 表明 SAVI 与 K 的 相关 性 最 好 。 本 文 在 
前 人 研究 的 基础 上 引入 SM 数据 ,通过 综合 作物 信 
息 .土壤 信息 ,使 建立 的 作物 系数 模型 更 加 可 靠 。 
3.3 K. 估 算 模 型 精度 验证 

为 评估 天 .估算 模型 的 可 行 性 ,将 2018 年 生长 季 
模拟 所 得 的 K. 值 乘 以 ET, 值得 到 E7, 的 模拟 值 ,以 涡 
度 相 关系 统 实测 实际 蒸 散 量 为 标准 进行 线性 回归 
分 析 , 由 于 草 旬 试验 区 2018 年 生长 季 内 5 月 涡 度 数 
据 缺 失 , 故 用 6 月 1 日 至 9 月 25 日 数据 进行 验证 ,并 
用 RMSE |d „bias 和 Nash-Sutcliffe 效 率 指数 4 个 统计 
量 对 模型 的 精度 进行 进一步 评价 。 

由 图 4 可 以 看 出 , 草 旬 和 沙丘 验 区 的 ET 模拟 值 


表 4 K. 估 算 模 型 


Tab.4 Linear equations for estimating corp coefficient(K.) 


站 点 变量 方程 RMSE Raj 
" NDVI(casel ) K,-1.287NDVI40.307 0.12 0.52 
= NDVI,SM(case2) K.=0.851 NDVI+3.108SM-0.640 0.06 0.84 
— LAI(casel ) K,-0.315LAI«0.098 0.14 0.57 

LAI,SM(case2) K.=0.239LAI+5.295SM-0.034 0.12 0.71 
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图 4 草包 和 沙丘 试验 区 模拟 蒸 散 量 和 实测 燕 散 量 季节 变化 及 回归 分 析 


Fig. 4 Seasonal variation and scatter plots of observed and estimated actual evapotranspiration at the meadow and sand dune locations 


与 实测 值 变化 趋势 基本 一 致 ,降雨 期 略 有 偏差 ,这 
是 由 于 雨 期 涡 度 相关 数据 精度 较 低 ,导致 模拟 的 
7 与 观测 到 的 有 7 之 间 存 在 差异 。 草 旬 试 验 区 生 
长 旺盛 季 蔡 散 强 烈 , 蔡 散 波动 较 大 , 涡 度 相关 测量 
时 段 内 的 实测 总 ET 为 372 mm, 模 拟 总 ET, 为 408 
mm, 日 平均 实际 蒸 散 量 为 3.38 mm:d', 模 拟 值 为 
3.66 mm .d ,高 估 了 8.28% , 偏差 较 小 ; 沙丘 试验 区 
生长 季 内 莹 散 量 大 致 呈 抛 物 线 型 曲线 变化 ,测度 实 
测 总 蒸 散 量 为 236 mm, 模 拟 值 为 261 mm ,高 估 了 
10.59% ;由 模拟 蒸 散 量 和 实测 蒸 散 量 回归 分 析 可 以 
看 出 , 草 负 和 沙丘 试验 区 的 模拟 值 与 实测 值 的 分 布 
均 较 为 集中 ,大 多 分 布 在 1:1 直线 附近 ,决定 系数 户 
分 别 为 0.76 和 0.57 ,模拟 效果 较 好 。 

由 表 5 可 知 ,RMSE 可 衡量 模拟 值 与 实测 值 之 间 
的 偏差 , 草 甸 和 沙丘 试验 区 的 RMSE 分 别 为 0.59 和 
0.42 ,表明 模拟 值 与 实测 值 较 接近 ;d 值 分 别 为 0.92 
和 0.85, 均 趋 近 于 1, 表 明 模 型 模拟 效果 好 ,bias 反 映 
实测 值 与 模拟 值 的 平均 偏差 ,两 试验 区 bias fH 
(-0.39、-0.11) 都 接近 0, 精 度 较 高 ,而 了 (0.64 和 
0.56) 均 接近 1, 表示 模型 可 信 度 高 ,综合 表 5 可 知 ， 
草 甸 和 沙丘 试验 区 均 有 较 好 的 模拟 效果 。 


表 5 作物 系数 到 估算 模型 的 适应 性 分 析 
Tab.S Adaptability analysis of K. estimation model 


评价 指标 RMSE d bias E 
px 0.59 0.92 -0.39 0.64 
沙丘 0.42 0.85 -0.11 0.56 


本 文 利 用 卫星 乏 感 数据 获取 植被 信息 ,结合 地 
面 实测 气象 .土壤 数据 建立 了 天 佑 算 模型 ,并 通过 
涡 度数 据 证 实 K. 估 算 模 型 在 干旱 区 草 甸 及 沙丘 均 
有 较 好 的 适用 性 。 但 本 文 并 未 考虑 植被 不 同 生育 
期 作物 系数 变化 情况 ,在 后 续 的 研究 中 将 结合 多 年 
数据 ,进一步 优化 到 估算 模型 ,深入 探讨 具体 生长 
阶段 内 作物 系数 的 变化 规律 及 其 影响 因素 。 


4 结论 

AS X fi] Hj Landsat 8 卫星 遥感 数据 ,人 研究 了 利用 
归 一 化 植被 指数 (NDVI) .土壤 调整 植被 指数 (SA- 
V 了 DD) 和 简单 比值 植被 指数 (SR) ,结合 地 面 实测 叶 面 
积 指数 (LAIT) 及 表层 土壤 含水 量 (SM ) 估 算 科 和 尔 沁 区 
草 旬 及 沙丘 试验 区 作物 系数 玉 的 可 行 性 ,得 到 以 下 
结论 : 

(1) 对 比 2017 年 生长 季 作 物 系数 K.、 植 被 指数 
及 SM 的 变化 趋势 发 现 ,各 植被 指数 和 SM 的 变化 规 
律 与 作物 系数 到 有 很 高 的 一 致 性 ,确立 基于 植被 指 
数 和 SM 估算 作物 系数 到 具有 一 定 的 可 行 性 。 

(2) 相关 性 分 析 崭 除 相 关 性 较 差 的 变量 SR, 逐 
步 回 归 分 析 再 次 剔除 不 显著 变量 ,最终 结合 土壤 水 
分 信息 de HEM LAIL X Landsat 8 遥感 数据 计算 的 
NDVI 建 立 玉 估算 模型 ,经 验证 ,模型 的 可 信和 度 较 
高 , 草 甸 和 沙丘 试验 区 的 尼 分 别 为 0.84 和 0.71。 

(3) 利用 2018 年 数据 对 作物 系数 到 估算 模型 
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进行 进一步 验证 , 草 铅 和 沙丘 试验 区 蒸 散 量 的 模拟 
值 与 实测 值 均 分 布 在 1:1 直线 附近 , 尼 较 高 (0.79 和 
0.56) ,统计 结果 表明 模拟 值 与 实测 值 有 较 好 的 对 应 
关系 [*] 
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Crop coefficient simulation and evapotranspiration estimation of dune 


and meadow in a semiarid area 


LI Xia, LIU Ting-xi", DUAN Li-min", TONG Xin", WANG Guan-li 
(1. Water Conservancy and Civil Engineering College, Inner Mongolia Agricultural University, Hohhot 010018, 
Inner Mongolia, China; 2. Inner Mongolia Key Laboratory of Water Resource Protection and Utilization, Hohhot 
010018, Inner Mongolia, China) 


Abstract: The crop coefficient is of great significance for improving the estimation accuracy of the actual evapo- 
transpiration and the regulation of regional water resources. In this paper, based on the Landsat 8 data extraction,we 
use three vegetation indices, normalized difference vegetation index (NDVI), soil adjusted vegetation index (SAVI), 
and simple ratio vegetation index (SR), combined with the ground measured soil moisture (SM) and the leaf area in- 
dex (LAI) data in 2017.We then use the stepwise regression analysis to establish a crop coefficient K. estimation 
model for the meadow and dune experimental area in the Horqin area. The simulated evapotranspiration(E7,) is ob- 
tained by multiplying the simulated K. value and the potential evapotranspiration (ETo) calculated by the FAO 56 
PM model in 2018. The accuracy of the K. estimation model was verified by comparing it with the actual evapotrans- 
piration(ET.) measured by the eddy covariance system. The results showed that the variation trend of K. in the mead- 
ow and dune experimental area during the growing season was consistent with VIs and SM, indicating that it was fea- 
sible to establish a K. estimation model based on these indicators. In correlation analysis, the correlation between K. 
and SR in the meadow experimental area was not significant, whereas the correlation between K. and SR in the dune 
experimental area was low. In the stepwise regression analysis, the insignificant factors were further removed, and 
the crop coefficient estimation model of meadow and dune was established; Ris was 0.84 and 0.71. The ET, calculat- 
ed on the basis of the K. estimation model achieved good simulation results both in the meadow and dune experimen- 
tal areas. 


Keywords: Horqin sandy land; evapotranspiration; remote sensing; crop coefficient; vegetation index; vorticity 


